arXiv:1505.03346v2 [cs.IT] 29 Dec 2016 


Analytic Solutions to a Marcum Q—Function-Based 
Integral and Application in Energy Detection of 
Unknown Signals over Multipath Fading Channels 


Paschalis C. Sofotasios^’^, Mikko Valkama^, Theodores A. Tsiftsis^’^, Yury A. Brychkov^, 

Steven Freear^ and George K. Karagiannidis^’® 

^Department of Electronics and Communications Engineering, Tampere University of Technology, 33101 Tampere, Einland. 

e-mail: {paschalis.sofotasios; mikko.e.valkama} @tut.fi 

^Department of Electrical and Computer Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece. 

e-mail: geokarag@ieee.org 

^Department of ElecUical Engineering, Technological Educational Institute of CenUal Greece, 35100 Lamia, Greece. 

e-mail: tsiftsis@teilam.gr 

^The Dorodnicyn Computing Center of the Russian Academy of Sciences, 119333 Moscow, Russia. 

e-mail: brychkov@ccas.ru 

® School of Electronic and Electrical Engineering, University of Leeds, LS2 9JT, Leeds, United Kingdom. 

e-mail: s.freear@leeds.ac.uk 

^Department of Electrical and Computer Engineering, Khalifa University, PO Box 127788, Abu Dhabi, UAE. 


Abstract —This work presents analytic solutions for a useful 
integral In wireless communications, which Involves the Marcum 
Q—function In comhination with an exponential function and 
arbitrary power terms. The derived expressions have a rather 
simple algebraic representation which renders them convenient 
both analytically and computationally. Furthermore, they can 
be useful in wireless communications and particularly in the 
context of cognitive radio communications and radar systems, 
where this Integral is often encountered. To this end, we derive 
novel expressions for the probability of detection in energy 
detection based spectrum sensing over rj—fj, fading channels. 
These expressions are given in closed-form and are subsequently 
employed in analyzing the effects of generalised multipath fading 
conditions in cognitive radio systems. As expected, it is shown 
that the detector is highly dependent upon the severity of fading 
conditions as even slight variation of the fading parameters affect 
the corresponding performance significantly. 

I. Introduction 

The generalized Marcum Q—function, Qm{a,b), is a vital 
special function in wireless communication theory. It was 
proposed several decades ago and has appeared extensively 
in various analyses in the context of stochastic processes in 
probability theory, single- and multi-channel based commu¬ 
nications over fading channels, information-theoretic analysis 
of multiple-input-multple-output (MIMO) systems, cognitive 
radio and radar systems, among others Q-®, and the refer¬ 
ences therein. Its use has also led to the derivation of numerous 
tractable analytic expressions, while its computational realiza¬ 
tion is rather straightforward since it is included as a built-in 
function in the most popular software packages inoi-iisi. 


However, it is widely known that the derivation of tractable 
analytic expressions in natural sciences and engineering can 
be rather laborious and cumbersome, if not impossible, as 
integrals that involve combinations of elementary and special 
functions are often required to be evaluated analytically oa- 
Ea, and the references therein. This is also the case when 
the Marcum Q—function is involved in integrands along with 
exponential and arbitrary power terms. A general form of such 
an integral is the following: 

poo 

Ia,b{k,m,p)= / x'^''~^Qm{ax,b)e~P'^^dx ( 1 ) 

Jo 

which can be equivalently expressed as 
1 

Ia,b{k,m,p) = - x’‘~^Qm{as/x,b)e~^''dx. ( 2 ) 

^ Jo 

The integrals in o and (|2]) are encountered in various 
applications relating to wireless communications, such as 
in the analysis of multichannel diversity systems with non¬ 
coherent and differentially coherent detection and in sensing 
of unknown signals in the context of cognitive radio and radar 
systems ll22l - ll40l and the references therein. Based on this, 
a recursive formula restricted to only integer values of k and 
m was reported in 0 while an infinite series representation 
for the case that k is arbitrary and m is positive integer was 
recently reported in IHTI . 

Nevertheless, these expressions are neither generic, nor ac¬ 
count for the case that m is an arbitrary real. Motivated by this, 
this work is devoted to the derivation of analytic expressions 


for Ia^b{k,m,p) which can be useful in applications relating 
to the wide field of digital communications. To this end, novel 
analytic expressions are derived for the probability of detection 
of energy detection based spectrum sensing over generalized 
rj—p fading channels. The derived expressions are given in 
closed-form and are utilized in analyzing the performance of 
the detector in various fading conditions. 

The remainder of this paper is organized as follows: Sec. 
II provides the derivation of two analytic expressions for 
Ia.bik,m,p). Sec. Ill is devoted to the application of the 
offered results in the analytic performance evaluation of energy 
detection based spectrum sensing over rj—p fading channels 
for various severity scenarios. The corresponding numerical 
results are given in Sec. IV along with useful discussions, 
while closing remarks are provided in Sec. V. 


II. Analytic Solutions to la, b{k,m,p) Integrals 

A. A Closed-Fonn Expression for Arbitrary Integer Values of 
k and Arbitrary Real Values of m 

As already mentioned, no analytic expressions for O and 
(|2|i for the case of arbitrary real values of m have been reported 
in the open scientific and technical literature. 


Theorem 1. For a,b,m,p G K'*' and k G N, the following 
closed-form expressions is valid. 


Ia,b{k,m,p) 


r(fc)r (m, 
2p^T{m) 


k-l 


+E 

1=0 




(3) 


where r(a:) and r(a;, a) denote the gamma function and upper 
incomplete gamma function, respectively, x\ = r(a; — 1 ) is the 
increasing factorial and iFi{x, y, z) is the Kummer confluent 
hypergeometric function h42’^ - f4^ . 


Proof: By integrating (|2|i by parts one obtains, 
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Xa,b{k,m,p) = lim 


f k—1 — 


P^dx 


J 0 


(4) 
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dx 


Qm{ay/x, b)dx 


where c is a non-negative finite real. By recalling that the lower 
incomplete gamma function is given by y{a,x) = r(a) — 
r(a, t), it readily follows that J x°‘~^exp{—x)dx = 7 ( 0 , x) = 
—r(a,x). Upon substituting in (01 one obtains. 


T = lim 

C—>-oo 


Qmiay/x, b)y{k,px) 


2p'^ 


J 0 


Qmi0,b)T{k,0) 

2p^ 


— lim 

c—>-oo 


Qm{as/c, b)T{k,pc) 

2p^ 


(5) 


With the aid of the identities for Qm{a,b) and r(a,x) func¬ 
tions in [?], IJI and ll42l as well as expressing dQra{a.,h) / da 


according to m eq. (10)], it follows that 


m—1 




l\p^2^ 




1=0 
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dy. 


By making the necessary variable transformation in 03] II. 7 
- pp. 726] and substituting in (|6|l yields 


Ia,b{k,m,p) = - 

1=0 

^ b^T{k)e-^ r 


(7) 


Notably, the above integral can be expressed in closed-form 
with the aid of 03l eq. (2.15.5.4)]. To this effect, by perform¬ 
ing the necessary change of variables and substituting in (|7]l, 
equation (O is deduced, which completes the proof. ■ 


It is noted that the algebraic representation of the derived 
solution is simpler than the recursive expression in 0, 11 
and additionally, the value of m is subject to no restrictions. 


B. An Exact Infinite Series toXa,b{k,m,p) for Arbitrary Reals 

It is recalled that no analytic expressions exist for dD and 
(|2]l for arbitrary real values, i.e. unrestricted, of all involved 
parameters. 


Lemma 1. For a, b, k,m,p G the following exact infinite 

series representation is valid for the integral in dB and d3. 


Xa,b{k,m,p) = ^ 
1=0 


a^^2'^T{k + l)T (m + l,^') 
l\T{m + 1) (a2 -h 


( 8 ) 


Proof: The Qm{a, b) function can be expressed in infinite 
series according to 16] eq. (29)]. Therefore, by performing the 
necessary change of variables it immediately follows that, 

Q„(aVS.b) = e-^-£ + W 

which upon substitution in dll yields. 


Xa,b(k,m,p) = 


00 (^m + l,^^ 

E l!2ir(m + l) 


00 

^(p+4) 


dx ■ 


n 


( 10 ) 


The above integral can be expressed in terms of the r(.) 
function in l43l eq. (2.10.3.2)], yielding 


2'=+'r(fc-f 0 
-I- 2p 


( 11 ) 


Evidently, by substituting (fTTl i in (ITOl i. equation d8]l is deduced 
thus completing the proof. ■ 





























The series in ([8]l is convergent and can provide acceptable 
accuracy when truncated after relatively few terms. However, 
deriving a closed-form expression for the truncation error is 
particularly advantageous in determining the corresponding 
truncation error accurately and straightforwardly. 


Lemma 2. For a,b,k,m,p G K'*", the following inequality 
can serve as a closed-form upper bound for the truncation 
error of 

_^r(fc) " a2'2'=r(A: + 0r(m + /,^) 

^ ;!r(TO + Z) (a2 + 

Proof: The truncation error of ® when this is truncated 
after n terms is expressed as, 

°° a2'2'=r(/c + Z)r(TO + Z,^) 

Z!r(m + Z) {a^ + 2pf+^ 
a2'2'=r(fc + Z)r (to + Z,4) 


= E 


\k+l 


(13) 


Z!r(m + Z) (a^ + 2 p) 
a 2 ' 2 '=r(Zc + Z)r (to + Z,^) 


-E' 

/=0 


Z!r(m + Z) (a2 + 2p) 


k+l 


It is recalled that the r(a,a;) function is monotonically de¬ 
creasing w.r.t. X and thus, r(a,a;) < r(a). To this effect, the 
upper incomplete gamma function in (fOT l can be bounded as, 

a^'- 2 ^T{k Fl) A °'^' 2 '=r(Z(; -f Z)r (m -f Z, 


1=0 


Z! (a^ -f 2p) 


1=0 


l\T{m Z) (a^ -I- 2p) 


k+l 


(14) 

By recalling the Pochhammer symbol, (a)n = r(a + n)/r(a), 
it follows that r(k + Z) = (Zc);r(Zc). Based on this, the above 
infinite series can be expressed as follows. 


„2lr,k 
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a2‘2'=r(fc-f Z) 
\T+T 


r(Zc)2'= 
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a^‘(k)i 


i. / 7- (15) 

(a^ -P 2p) Z! (a^ + 2p) 
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The infinite series in the right-hand side of (HI can be 
expressed in terms of the hypergeometric function, namely. 


a^‘(k)i(l)i 
^ Z! (a2 -f 2p)‘ {l)i 
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= i-Z^o fc;; 


2p 


Based on (fTST l. it immediately follows that 


1^0 k] ; 


2p 


(a^ -I- 2pY 

2kpk 


Therefore, by substituting in ( fTTI i in (flSl l one obtains. 


“ a^^2'^T{k + l) T{k) 


E 


^ Z! (a2 -f 2p) 


k+l 


(16) 


(17) 


(18) 


Evidently, by substituting (fTsT i in (fT4li yields (fTSli thus, com¬ 
pleting the proof. ■ 

To the best of the Authors’ knowledge, equations ®, 
® and (fTST i have not been previously reported in the open 
technical literature. 


III. Applications in Cognitive Radio 


A. Energy Detection Based Spectrum Sensing 

Cognitive radio (CR) is an emerging technology that allows 
opportunistic access of licensed frequency bands when they 
are not utilized. Given the increased spectrum scarcity along 
with the high demands for bandwidth resources, CR is antic¬ 
ipated to play a core role in the next generation of mobile 
communication systems, namely 5G. The most important part 
of CR technology is the accurate and robust sensing of vacant 
frequency bands and based on the respective decision the user 
will decide on whether it can establish communication or 
not. Therefore, spectrum sensing is the most critical operation 
in CR systems with energy detection being regarded as the 
most simple and popular method ll2Tl . In this context, the 
performance of energy detection based spectrum sensing over 
various fading conditions have been investigated in ll24ll - ll^ . 
llTSll - and the reference therein. 

It is recalled that in narrowband energy detection, the 
received signal waveform follows a binary hypothesis that can 
be represented as ll^ eq. (1)], 


r{t) = 


n{t) 

hs{t) n{t) 


■■Ho 

:iJi 


(19) 


where s{t), h and n{t) denote an unknown deterministic 
signal, the amplitude of the channel coefficient and an additive 
white Gaussian noise (AWGN) process, respectively. The 
samples of n(t) are assumed to be zero-mean Gaussian random 
variables with variance NqW with W and Nq denoting the 
single-sided signal bandwidth and a single-sided noise power 
spectral density, respectively lf3^ . The hypotheses Hq and 
Hi refer to the cases that a signal is absent or present, re¬ 
spectively. The received signal is subject to filtering, squaring 
and integration over the time interval T which is expressed 
as || 22 l eq. ( 2 )], y = fo I P dt. The output of 

the integrator corresponds to a measure of the energy of the 
received waveform and acts as a test statistic that determines 
whether the received energy measure corresponds only to the 
energy of noise {Ho) or to the energy of both the unknown 
deterministic signal and noise {Hi). By denoting the time 
bandwidth product as u = TW, the test statistic typically 
follows the central chi-square distribution with 2u degrees of 
freedom under the Ho hypothesis and the non central chi- 
square distribution with 2u degrees of freedom under the Hi 
hypothesis ll24l . Based on this and by recalling that energy 
detection is largely affected by a predefined energy threshold, 
A, the performance of the detector is characterized by the 
probability of false alarm, Pf — Pr(j/ > A | Hq) and the 
probability of detection, Pd = Pr{y > A | Hi), namely ll22l . 


Pf = 


r(»71) 

r{u) 


( 20 ) 


and 


Pd = Qu 



( 21 ) 


respectively. 



















B. The ri—fj, Distribution 

The T]—^ distribution is a generalized fading model that 
has been widely shown to provide adequate characterization 
of multipath fading in non-line-of-sight (NLOS) communi¬ 
cations. It was reported in ||45]| along with the n—fi fading 
model which accounts for corresponding line-of-sight (LOS) 
communication scenarios. The 77 —p fading model has been 
shown to be particularly flexible and it includes as special 
cases the well known Hoyt, Nakagami— to, Rayleigh and one¬ 
sided Gaussian distributions Il45l . Its remarkable flexibility and 
usefulness were demonstrated clearly in m Fig. 9] along 
with the K—fi fading model where it is clearly shown that the 
77—77 fading model is significantly more flexible than the more 
commonly adopted Nakagami— to and Rayleigh distributions. 

In terms of physical interpretation, the rj—fi fading model 
is expressed by two physical parameters, 77 and p, and it 
holds for two formats, namely Fonnat-1 and Fonnat-2. In 
the former, the 77 parameter denotes the ratio of the powers 
between the multipath waves in the in-phase and quadrature 
components, whereas in the latter it denotes the correlation 
coefficient between the scattered wave in-phase and quadrature 
components of each cluster of multipath. Likewise, the p 
parameter denotes - in both formats - the inverse of the 
normalised variance and relates to the number of multipath 
clusters in the environment ||45]| . 

The SNR probability density function of the 77 —p distribu¬ 
tion is expressed as. 


^ 7 ( 7 ) 


r(p)iT'"-5 ^ J 


( 22 ) 


where 7 denotes the average SNR whereas 

2 + i]-^+r] r]-^-T] 

- 7 -) H = --- 

4 4 


(23) 


in Format-1 with 0 < rj < 00 and, 

h=-^, H = -^ 

1 — 77^ 1 — 77^ 

in Format-2 with — 1 < 77 < 1. In addition, 

L;2(i?2) r H 

^ ~ 2Var(i?2) [ X 


(24) 


(25) 


with E{.) and f denoting expectation and the root-mean- 
square (rms) value of the envelope R, respectively na. 


C. Energy Detection over rj—gi Fading Channels 

Corollary 1. For 77 , 7 , A € M'*" and p G N, the following 
closed-form expressions hold for the average probability of 

*The Format -2 of the 77—77 distribution is also known as A—77 distribution. 


detection over rj—p fading channels, 

f i-iy i-ir \ 

^ i\2F+iHF+i \{h- ny-i {hF ny-i j 

i —0 

(—1)* iFi ^1 -I- 7 , 1 -I- 77 , 27 + 4 / 777 - 77 ) ) 

{h - F 2{h - iF)p)*+i 

( — 1)^ iFi -f 7 , 1 -I- 77 , 27 + 4 /tp 7 + 77 )) 

{liF F2{hF H)py+^ 

(26) 


where g{a,x) = r(a,x)/r(a) denotes the regularized upper 
incomplete gamma function and h and H are given by (1231) 
and (El according to Format-1 and Format-2, respectively. 


Proof: The average detection probability is obtained by 
averaging (|2TI) over the fading statistics of the channel, namely, 

= y Qu v^) P-fh)dy. (27) 

By substituting (l22ll in dZTl l one obtains. 



where 

~ r(p)i7'^-^7'"+5 ■ 


(28) 

(29) 


Importantly, for the special case that p is a positive integer, 
the Bessel function in (|28] | can be expressed in closed-form 
with the aid of ll46l eq. (8.467)] namely, 

I , = V {-lyrjpFiW+h^-^ 

7 J t^(!v^r(p-()(4pi77)'+^ 

77 — 1 7 , 1 2uH-y 

^ (^^ (-l)Mr(p + 0f+ie-^ 

Uy/^T{p-l){ApHyy+^ 

Therefore, by substituting (1^ in (l28l l it follows that. 


p A{-iyr{pFiW+i 

^ (! 0 Fr(p-()( 4 piT )'+5 

^ A{-irripFiW+^ 
^(!0Fr(p-O(4pi7)'+^ 



Qu (a/27) a/a) dy 


yZ —/i + 1 


2M(h-tf)7 




Notably, the integrals in (OTT i have the same algebraic represen¬ 
tation as ([T]) and (|2| and thus, they can be expressed in closed- 
form with the aid of Theorem 1. As a result, by performing the 
necessary change of variables in (12), substituting in (|3T1) and 
carrying out long but basic algebraic manipulations, equation 
(l26l ) is deduced and thus, the proof is completed. ■ 
































Remark 1. The energy threshold in can be expressed 
as X = 2Q~^ {u,Pf), where is the inverse regularized 

upper incomplete gamma function. To this effect (I261 l can be 
also equivalently expressed in terms of Pf as follows: 


fi-i 




{T)ih^Pf 

\ {h-HY-’- ' {h + Hy-^ 


1=0 
/i— 1 ^s—l—1 

E E 

l—O i=0 


(- 1 )' 


(-1)^ 




(-1)' iFi (^1 + 1 + u, 

{h - + 2{h - H)p,y+^ 


{-ly iFi (^1 + i, 1 + M, ) 

{h + Hy-^-^iff + 2{h + 


(32) 


To the best of the Authors knowledge (l26l l and (l32l l have not 
been previously reported in the open technical literature. 


IV. Numerical Results 

This section is devoted to the analysis of the behaviour 
of energy detection in p—p. fading conditions by means of 
Pd versus 7 curves and complementary receiver operating 
characteristics (ROC) curves {Pm versus Pf). To this end, Fig. 
1 illustrates the behavior of the Pj. versus 7 for different values 
of the fading parameters p and p for constant time-bandwidth 
product u = 3 and the case that Pf = 0.01 and Pf = 0.1. One 
can notice that the average probability of detection is increased 
as p increases from 0.01 to 0.95 for both cases of Pf. This 
is also the case for the p parameter as for a fixed value of 
p, the Pd increases when p = 3 compared to the case that 
p = 1. This also holds for both Pf = 0.01 and Pf = 0.1 and 
particularly for moderate to high average SNR levels. 

In the same context. Fig. 2 depicts the corresponding ROC 
curves {Pm = 1 — Td versus Pf). The value of Pf is assumed 
between 0.01 and 0.2 while u = 4 and 7 = 15dB. One 
can observe how the performance of the detector improves 
as the severity of fading is reduced in terms of both p and p. 
Indicatively, for Pf = 0.1 the value of Pm reduces by over 
70% when p changes from 0.01 to 0.95 for p = 1.0 and 
over 65% when p changes from 1.0 to 2.0 for p = 0.95. This 
demonstrates the sensitivity of the energy detector in multipath 
fading conditions and how the corresponding severity of fading 
can affect its performance and robustness. 


V. Conclusion 

New expressions were derived for a Marcum Q—function 
based integral that is often encountered in the broad area of 
digital communications. These expressions include an exact 
closed-form expression and an infinite series representation 
along with a closed-form upper bound for the corresponding 
truncation error. The offered results have a relatively simple 
algebraic representation which render them analytically and 
numerically convenient while they can be useful in numerous 
applications in wireless communications. As an example, they 



Fig. 1. Pd vs 7 for different values of rj and with u = 3 and Pf = 0.01 
and Pf = 0.1. 



Fig. 2. Pm vs Pf ROC curve for for different values of rj and with 
u = 4 and 7 = 15 dB. 


were used in energy detection based spectrum sensing, in 
the context of cognitive radio and radar systems, deriving 
novel closed-form expressions for the average probability of 
detection of unknown signals over p—p fading channels. The 
derived expressions were subsequently employed in analyzing 
the effect of multipath fading on the spectrum sensing perfor¬ 
mance and it was shown that the overall performance of the 
detector is, as expected, largely affected by the value of the 
involved fading parameters, particularly for moderate and high 
SNR levels. 
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